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Abstract 

X 

■ A lattice QCD calculation of the proton matrix element of the flavor 

singlet axial-vector current is reported. Both the connected and 
disconnected contributions are calculated, for the latter employing the 
variant method of wall source without gauge fixing. From simulations in 
quenched QCD with the Wilson quark action on a 16^ X 20 lattice at 
(3 = 5.7 (the lattice spacing a ~ 0.14fm), we find 

AS = An Ad As = -h0.638(54) - 0.347(46) - 0.109(30) = -h0.18(10) 
with the disconnected contribution to An and Ad equal to —0.119(44), 
which is reasonably consistent with the experiment. 
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The flavor singlet axial- vector matrix element of proton has been widely 
discussed in recent years. The interest initially arose from the EMC data[l] for 
the spin-dependent proton structure function gi which, taken together with 
earlier SLAC data[2], apparently indicated that the fraction of proton spin 
carried by quarks has a small value AE — Au + Ad -|- As = 0.12(17) and that the 
strange quark contribution is unexpectedly large and negative As = — 0.19(6) [1]. 
New experiments have since been performed with proton[3, 4], deuteron[4, 5] and 
neutron[6] targets. Combined reanalyses of these data using 
qa^ Au-Ad^ 1.2573(28) and Au + Ad- 2 As = 0.601(38) [7], with the 
aid of three-loop perturbative QCD calculations for the structure functions[8], 
have recently been carried out [9, 10], reporting for the matrix elements 
renormalized at /i^, 

AE = Au + Ad + As 

f 0.83(3) - 0.43(3) - 0.10(3) = 0.31(7), //^ = 10GeV^[9] 

\ 0.832(15) - 0.425(15) -0.097(18) = 0.31(4), //^ = oo[10]. ^ ^ 

where the first analysis[9] employed the Bjorken sum rule to fix the value of the 
strong coupling constant ag while the world average for as was used in the 
second analysis [10]. 

In this article we report [11] on a quenched lattice QCD calculation of the 
proton matrix elements of axial- vector current ior u, d and s quarks including 
both connected and disconnected contributions. A serious technical obstacle in 
such a calculation has been a reliable estimate of the disconnected part, which we 
are now able to overcome with the variant of the method of wall sources [12]. An 
exploratory study employing a random Z{2) source was previously made in 
Ref. [13]. Majority of lattice QCD calculations to date, however, attempted to 
evaluate AE from the proton matrix element of topological charge 
density [14, 15, 16]. The use of quenched approximation is not valid in this 
approach due to the degeneracy of rj' and 7r[14]. It also turned out that data 
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generated in quenched QCD[14, 16] are too noisy to extract the matrix element. 
Prom a full QCD calculation with four flavors of Kogut-Susskind quarks, the 
authors of Ref. [15] reported AE = 0.18(2). This is a difficult calculation and the 
quahty of their raw data does not seem to be as good as the error they quoted 
indicates. A direct calculation of axial- vector current matrix elements is superior 
in that the problem due to ry' is absent and that the contribution of connected 
and disconnected contributions for each quark flavor and their quark mass 
dependence can be examined. 

Let us define < Ps|?7375?|Ps s • Aq with \ps > the proton state at rest with 
the spin projection in the z direction equal to s/2. To extract Aq we calculate 
the ratio of the three-point function of the proton and axial-vector current to the 
proton two-point function, each projected onto the zero momentum state, 

p/.N E.=±g- <Psit)Et'^0Ql3l5q{t')ps{0) > large t + , 7-I A . /o^ 

""^'^^ E...<P.(t)p.(0)> > const. + Z,Aqt, (2) 

with Za the lattice renormalization factor for the axial- vector current. The 
connected amplitude can be calculated by the conventional source method[17]. 
To handle the disconnected piece we employ quark propagators G{n,t) evaluated 
with unit source at every space-time site ( except for the t — time slice to avoid 
mixing with the proton valence quark propagator) without gauge fixing[12]. The 
product of the nucleon propagator and J2{n,tjto) T^[l3l5G{'n.,t)] equals the 
disconnected amplitude up to gauge- variant non-local terms which cancel out in 
the average over gauge configurations. 

Our calculation is carried out for the Wilson quark action in quenched QCD 
a,t P — 5.7 on a 16^ x 20 lattice. The u and d quarks are assumed to be 
degenerate with Ku = = Kq, while the strange quark is assigned a different 
hopping parameter Ks- We analyzed 260 configurations for the hopping 
parameters Kq, Kg — 0.160, 0.164 and 0.1665, generated with the single plaquette 
action separated by 1000 pseudo-heat bath sweeps. In order to avoid 
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contaminations from the negative-parity partner of the proton propagating 
backward in time, we employ the Dirichlet boundary condition in the temporal 
direction for quark propagators. We use the relativistic proton operator, fixing 
gauge configurations on the t = time slice to the Coulomb gauge in order to 
enhance proton signals. Let us emphasize that the gauge fixing is limited to the 
t — time slice; the rest of the lattice is left unfixed. Thus the issue of 
contribution of gluons to proton spin through gauge non-invariant Chern-Simons 
current [18] does not arise in our case. 

In Table 1 we list the results for hadron masses obtained by a single 
exponential fit of propagators over 6 < t < 12. Errors in this table and below are 
estimated by the single elimination jackknife procedure. 

We plot the connected contribution of u and d quark to the ratio R{t) in 
Fig. 1 (a) and (b) for the case of Kg = 0.164. Very clean signals with the linear 
behavior up to i w 14 are observed, with the sign of the slope consistent with 
expectations from quark models which predict Au = 4/3 and Ad — —1/3 in the 
static limit. For the disconnected contribution the quality of our data is not 
quite good, in spite of reasonably high statistics of the simulation ( Fig. 1 (c) ). 
A region showing a linear dependence is limited to i ~ 5 — 10 and errors grow 
rapidly with increasing t; For t > 12 signals are lost into a large noise. 
Nonetheless, we can still observe the negative value of the sea quark contribution 
to the axial- vector matrix element including that of strange quark, which is not 
predictable in quark models. To extract the axial-vector matrix elements, we fit 
the data for R{t) to the linear form (2) with the fitting range chosen to be 
5 <t < 10. Changing the fitting range to5<i<llor6<i<10 increases the 
value of fitted slope by 30 — 40% with a roughly proportional increase of error. 
(Including larger values of t does not seem reasonable due to a departure from a 
linear behavior for i > 11 and a rapid loss of signals.) 

Results are corrected by the tadpole-improved renormalization factor given 
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by[19] 

Za^{1- 0.31a,) (l - , (3) 

where we use a-jjgil/a) = 0.2207 for a,. We should note that the flavor singlet 
axial-vector current requires an additional lattice-to-continuum divergent 
renormalization from diagrams containing the triangle anomaly diagram. We 
leave out this factor since the explicit form of this contribution which starts at 
two-loop order has not been computed yet. The results for Au, Ad and As are 
tabulated in Table 2. 

We present the axial-vector matrix element for u, d and s quarks in Fig. 2 as 
functions of the bare u and d quark mass rUqa — {1/ Kq — l/Kc)/2 using 
Kc — 0.1694. For the strange quark the values represent result of an 
interpolation to the physical strange quark mass corresponding to m/^/mp = 0.64 
as explained below. As we already remarked the disconnected contributions 
(filled symbols) are negative. Furthermore their magnitude is small, albeit 
increasing slightly toward the chiral limit. The connected contributions exhibit 
an opposite trend of a slight decrease. 

We calculate the physical values of matrix elements in the following way. For 
u and d quarks, wc estimate the sum of disconnected and connected 
contributions by first combining the two contributions in the ratio R{t) and then 
fitting the result to the hnear form (2) over 5 < t < 10 for each Kq. The fitted 
values are extrapolated linearly to the chiral limit rUq = 0. For the strange quark 
contribution similar extrapolations to rUq — are made for each Kg, whose 
results in turn are fitted to a linear function of the strange quark mass 
TTiga = {l/Ks — 1/ Kc)/2. The physical strange quark mass is estimated by 
generalizing the relation m^a? — 2.710(25)mqa, obtained from the hadron mass 
results in Table 1, to m^a? — 2.710(25)(mga -|- msa)/2 and using the 
experimental ratio iriK/Tnp = 0.64, which gives m^a = 0.0826(22) or = 0.1648. 
This analysis, as summarized in Table 3, yields for the quark contribution to 
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proton spin, 

AS = Am + Ad + As = +0.638(54) - 0.347(46) - 0.109(30) = +0.18(10), (4) 

These values, notably the sign and magnitude of the strange quark contribution, 
show a reasonable agreement with the phenomenological estimate (1)[10]. 

Table 3 also hsts the flavor non-singlet matrix elements Fa — (Ait — As)/ 2 
and Da = (Am — 2Ad + As)/2 where the strange quark mass is interpolated to 
the physical value. In the chiral limit for u and d quarks niq = we find 
Fa = 0.382(18) and Da = 0.607(14), which implies qa^Fa + Da^ 0.985(25) 
and Fa/Da = 0.629(33). Compared to the experimental values qa = 1.2573(28) 
and Fa/ Da — 0.586(19) [7], the ratio shows a good agreement while the 
magnitude of Fa and Da are about 25% smaller. Previous quenched results at 
13 = 6.0[20, 21] and with the VS-blocked Wilson action[22], and for fuU QCD[20] 
at P — 5.4 — 5.6, with the lattice spacing of a ~ 0.15 — O.lfm, yield similar results 
if analyzed with the same renormalization factor as we employed. 

Possible sources of systematic errors in our results are scaling violation effects 
due to a fairly large lattice spacing a ~ 0.14fm of our simulation at P — 5.7 and 
uncertainties in the perturbative estimate of the renormalization factor (3). The 
small values of flavor non-singlet couplings compared to experiment by about 
25%, possibly arising from these uncertainties, suggest that our result for AE 
might be underestimating the continuum value by a similar magnitude. The use 
of quenched QCD might also cause systematic errors. We note that the lack of 
two-loop calculation for the flavor singlet lattice-to-continuum renormalization 
factor makes it difficult to specify the scale at which AE is evaluated, although 
we expect the scale dependence to be weak being a two-loop effect. While these 
points should be examined in future studies, we feel that the encouraging result 
we found points toward an eventual resolution of the spin crisis issue within 
lattice QCD. 
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Let us finally note that a result similar to ours has recently been reported at 
the Lattice '94 Symposium [2 3]. Employing the Z{2) noise method for evaluating 
the disconnected contribution on 24 configurations for an 16^ x 24 lattice at 
P = 6.0 in quenched QCD, the authors found 
AE = Axi + Ad + As = +0.78(7) - 0.42(7) - 0.13(6) = +0.22(9). 
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Figure Captions 

Fig. 1 Connected and disconnected contribution to R{t) for u and d quarks at 
P = 5.7 and Kg = 0.164 on an 16=^ x 20 lattice. 

Fig. 2 Axial-vector matrix elements for u, d and s quark as a function of 

degenerate u and d quark mass. Strange quark mass is interpolated to the 
physical value. 

Tables 



Table 1: Hadron mass results on an 16^ x 20 lattice at /3 = 5.7 obtained with 260 
configurations. 







rUpa 




0.1600 
0.1640 
0.1665 


0.6873(24) 
0.5080(29) 
0.3674(39) 


0.8021(29) 
0.6822(38) 
0.6085(58) 


1.2900(60) 
1.0738(80) 
0.915(11) 



Table 2: Nucleon axial- vector matrix elements as a function of K, 





Altconn. 


Adconn. 


Altdisc. = Addisc. 




As 












= 0.160 


0.164 


0.1665 


0.1600 


0.9071(92) 


-0.2470(35) 


-0.025(10) 


-0.025(10) 


-0.035(12) 


-0.042(15) 


0.1640 


0.839(19) 


-0.2382(87) 


-0.066(23) 


-0.049(19) 


-0.066(23) 


-0.076(27) 


0.1665 


0.818(39) 


-0.231(23) 


-0.093(54) 


-0.068(35) 


-0.084(46) 


-0.093(54) 


A", = 0.1691 


0.763(35) 


-0.226(17) 


-0.119(11) 


-0.083(31) 


-0.105(12) 


-0.117(50) 



Table 3: Total contribution to Am and Ad, and As at the physical strange quark 
mass corresponding to Kg = 0.1648. Flavor non-singlet coupling and are 
also listed. 





AWtotal 


Afitotal 


As 


Fa 


Da 


0.1600 
0.1640 
0.1665 


0.882(14) 
0.773(28) 
0.715(64) 


-0.273(11) 
-0.305(24) 
-0.326(56) 


-0.0374(89) 
-0.069(16) 
-0.087(33) 


0.4536(46) 
0.4196(93) 
0.409(20) 


0.7011(37) 
0.6580(71) 
0.641(16) 


= 0.1694 


0.638(54) 


-0.347(46) 


-0.109(30) 


0.382(18) 


0.607(14) 
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